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ABSTRACT 

We study single input and single output stochastic linear time-varying systems and develop a unified approach 
for generalized minimum variance controllers. The plant to be controlled is a time-varying controlled autoregressive 
moving average model and the performance index includes both tracking error variance of filtered plant output and a 
quadratic function of filtered plant input. Both the input and output filters are linear time-varying. 
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INTRODUCTION 

Generalized minimum variance controllers (GMVCs) are flexible controllers for optimal control of 
stochastic plants. They have seen many applications [1-3]. Various GMVCs have been developed based on different 
control strategies in order to meet various requirements of industrial applications [4-8]. The performance index of a 
GMVC is a cost functional that is a tracking error variance of filtered plant output penalized by a quadratic function of 
filtered plant input. The basic difference in different approaches is in the choice of input/output filters and weightings 
for the tracking error variance and the quadratic input function. The input and output filters can be used to modify the 
closed-loop dynamics for desired transient performance and the weightings can be used for specification of 
compromise between the output tracking accuracy and magnitude of fluctuation in plant input. 

In this paper, we develop a unified approach to various GMVCs for both linear time invariant (LTI) and 
linear time varying (LTV) single input and single output (SISO) systems by introducing general form of LTV filters 
for both plant input and output in the cost functional such that various GMVCs are special cases of this LTV GMVC. 
We will also provide closed-loop dynamics and stability analysis for the unified approach of GMVC. The plants to be 
controlled is described using an LTV controlled autoregressive moving average (CARMA) model and the filters are 
described using LTV transfer operators. 

The reminder of this paper is organized as the following. Section 2 describes the LTV SISO plants and the 
GMVC performance index. Section 3 develops the unified LTV GMVC and Section 4 presents a simulation example. 
Linally, Section 5 concludes this paper. 
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LTV PLANT AND CONTROL OBJECTIVE 

The LTV systems are described using the following SISO CARMA model. 

y{k+d) = A -1 (k, q~ l )[B{k, q~ l )u(k) + C(k , q~ l )w(k+d)] (1) 

where u(k) and y(kj are input and output, d is a positive integer representing the input/output delay, w(k) is an 
independent Gaussian noise with zero mean and a uniformly bounded time-varying variance, q 1 is the one-step-delay operator 
satisfying q' I f(k)g(k)=f(k-l)q' 1 g(k)=f(k-l)g(k-l). In the CARMA model A(k,q~ l ), B(k,q~ l ) and C(k,q -1 ) are LTV 
moving average operators (MAO's) having the following form. 

K(k,q- l ) = k 0 (k) + k l (k)q~ l +... + k nk (k)q~" k (2) 

where K(k,q~ l ) = A(k,q~'), B(k,q~ l ) and C(k,q~ l ) are polynomials in the delay operator q' 1 and 
a Jk) = c 0 (k) = 1 ■ The LTV MAO’s are natural extension of the moving average operators from the LTI transfer functions for 
LTV CARMA models, where the constant coefficients in the LTI polynomials are replaced using time-varying parameters. It 
is assumed that all the time-varying parameters of the LTV polynomials are uniformly bounded, i.e. \k (£)l<c<0 for a 

constant c and all the i and k. It is also assumed that the time delay between the plant input and output is time-invariant, i.e. 

I b 0 ( k ) l> c > 0, for a constant c and all k. The inversion operation of the MAO is an autoregressive operator (ARO) denoted 

by A^\k,q^ 1 ). A LTV ARO is also a natural extension of an LTI ARO of LTI transfer functions. Detailed properties of the 
LTV ARO can be found in [9]. For the design of the LTV GMVC it is assumed that both A^'(k.q^ 1 ) and C~\k,q~ ‘)are 
exponentially stable. 

Given a uniformly bounded reference z(k) the generalised minimum variance control objective is to find an LTV 
controller that generates a sequence of control variable u(k) such that the following generalised minimum variance cost 
functional is minimized. 

J(k+d) = E(\ P(k,q- 1 )y(k+d)-Q(k,q-')z(k) I 2 (3) 

+ V(k)\R(k,q-')u(k)\ 2 /D(k)) 

where D(k)={y(k), u(k), y(k-l), u(k-l), ...} is the set of input and output data up to and including the current time k 
and E(X/D(k )) is the operator for mathematical expectation of X conditioned on D(k), V(k)>c>0, for a constant c and all k, is 
for the weighting between the tracking error and the magnitude of input. The LTV MAO’s/^A:,^ 1 ), Q(k, q 1 ) and R(k,q ') 
are the LTV filters for the output, reference and input. They have the following forms. 

K(k,q~ l ) = K n (k,q~ 1 )K~\k,q~ 1 ) (4) 

where 

K{k,q-') = P(k,q ~ l ), Q{k,q~ l ), R(k,q~ l ) (5) 

and 
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K n (k,q~ l ) = k n0 (. k ) + k nl ( k)q ~ 1 (g) 

+t, 2 (t) 9 - j +...+fc w (t)r 

It is assumed that k n0 (k) is uniformly bounded away from zero and all the parameters in ( 6 ) and (7) are uniformly 

bounded away from infinite. The three filters are LTV moving average and autoregressive operators that are natural 
extensions of the LTI filters from LTI GMVCs for LTV GMVCs. 

Generalized Minimum Variance Controller 

We first develop a minimum variance predictor (MVP) for the filtered output 

y/{k + d) = P(k,q~ l )y(k + d). ( 8 ) 

The cost functional is 


J p (k + d) = E(\P(k,q~ 1 )y(k + d)-ys(k + d)\ 2 / D(k)) (9) 

Where y/(k+d) is the cZ-step-ahead prediction of if/(k). Left multiplying (1) using P(k,q ~ ') and noting ( 8 ) we 

have 


iff(k + d) = P(k,q ')A '( k,q l )B(k,q l )u(k) ^ 0 ^ 

+P n (k,q~ 1 )P; 1 (k,q- 1 )A-\k,q- 1 )C(k,q- l )w(k + d). 

The filtered output has two components. The first is the response to the input and the second is caused by the noise. 
The last term in the above equation is the component caused by the noise. 

Noting the last term in (10) and applying long division we have 

Pf\k,q-')A~\k,q~ l )C{k,q-')w(k + d) 

= [A(k,q~ 1 )P d (k,q~ 1 )Y 1 C(k,q~ 1 )w(k + d) ( 11 ) 

= {F(k,q- 1 ) + [A(k,q- 1 )P d (k,q- l )Y 1 G(k,q- l )q- d }w(k + d) 

= F(k,q- 1 )w(k+d) + [A(k,q~ 1 )P d (k,q- 1 )Y l G(k,q~ l )w(k) 

Where 

F(k,q~ l ) = l + f { (k)q- 1 + f 2 (k)q~ 2 +...+ f d _ t (k)q~ d+l (12) 

Is the quotient of the following long division 

[A(k,q~ l )P d (k,q~ l )]~ 1 C(k,q~ 1 ) 

= F(k,q- l ) + [A(k,q- 1 )P d (k,q- 1 )T l G(k,q- 1 )q- d 

Where 
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G(k,q l ) = g 0 (k) + g l (k)q~ 1 + ...+ g s (k)q~ 


(14) 


Is the remainder. The noise component can be further divided into two parts. The first is caused by the noise to occur 
in the future and the second is caused by the current and past noise. The future part is completely unknown and the current and 
past part can be estimated using the data set D(k). Left multiplying p (k.q ' ) on both sides of equation (1 1) we have 


P n ( k,q )P;' ( k,q _1 )A~ l (k , q~ l )C(k, (f 1 )w(k + cl) 

= P tI (k,q~ l )F(k,q~ 1 )w{k + d) 

+ P n (k,q- l )[A(k,q- 1 )P d (k,q- l )r 1 G(k,q- 1 )w(k). 


Letting 


P n (k,q~ 1 )F(k,q~ 1 ) = H(k,q~ 1 ) + L(k,q~ 1 )q~ d 


Where 

F[{k,q- x ) = l + h\{k)q~ x +... + h d _ l (k)q~ 


(15) 


(16) 


(17) 


Has all the terms that have delays less than d in P n (k,q ')F(k,q ') and L{k,q 1 )q d has all the rest terms, which 
have delays greater than d-1. Substituting (16) into (15) and noting (10) we have 


y/(k + d) = P(k,q ‘)H \k,q l )B(k,q l )u(k) 

+ H(k,q~ l )w(k + d) + \L(k,q~') 

+ P„ (k,q~ l )P;' (k, q- l )A~ l (k, q l )G(k , q~ l )]w(k) 


(18) 


Where H(k,q l )w{k + d) is the response due to future noise and [L(k,q 1 ) + P(k,q ‘)A 1 (k,q l )G(k,q 1 )\wik) 
has all the past and current noises for the filtered output. Taking mathematical expectation conditioned on D(k) on both sides 
of the above equation we have the following tZ-step-ahead minimum variance prediction of the filtered plant output. 


y/{k + d) = P(k,q *)A '(k,q l )B(k,q ')u(k) 

+ [L{k,q~') + P{k,q~ i )A~\k,q~ i )G{k,q~ l )]w{k) 
= P(k,q~ 1 )[A-\k,q~')B(k,q~ 1 )u(k) 

+ A’ 1 (k, q~ l )G(k, q^)w(k)] + L(k, q~ l )w{k) 

= P(k,ci l ) y(k + d) + L(k,q~ l )w(k) 

Where 


(19) 


y(k + d) = A l (k,q l )B(k,q l )u(k) 
+ A~ x {k,q~ l )G{k,q~ l )w{k) 


( 20 ) 


Is the t/-step-ahead minimum variance prediction of the plant output y(k) [10]. The above equation shows that the 
minimum variance prediction of the filtered output is not the filtered minimum variance prediction of the plant output. The 
difference between the two is L(k,q~ l )w(k) ■ 


Impact Factor (JCC): 5.7294 


NAAS Rating: 3.11 



A Unified Approach for Generalized Minimum Variance Controller for Linear Time-Varying Systems 


65 


GMVC Theorem 

If the LTV ARO’s A~ l (k,q~ l ), C~ l (k,q ~ l )& nd P~ l {k,q~ l ) are exponentially stable, the LTV GMVC for the 
CARMA model (1) is given by 


w(k) = C \k-d,q l )A(k-d,q l )y(k ) 

— C~\k — d,q~ l )B(k — d,q~ l )u(k) 

[B(k, q-' ) + A(k, q~ l )P~ l (k, q ~ l ) R(k, q~ l )]u(k) 

Pjk)b 0 (k) 

= A(k,q- l )p- l (k,q- l )Q(k,q- l )z(k ) 

-[A(k, q- 1 )P (k, q-' )L(k,q~ l ) + G(k , q~ l )]w(k). 

Proof. Subtracting (18) from (19) we have 

y/(k + d) = y/(k + d) - H(k,q~ l )W(k + d). 

Substituting (23) into the performance index (3) we have 

A 

J{k + d) =1 y/(k+d)-Q(k,q~ l )z(k) I 2 
+V(k) I R(k,q~ l )u(k ) I 2 +£(l H(k,q~ l )w(k + d) I 2 ). 

It follows that 

dJ ^ d) = 2 p n0 (k)b 0 mkk + d) 

ou(k) 

- Q(k,q- l )z(k)] + 2 r n0 (k)V(k)R(k, q~ l )u{k) 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 


and 


d 2 J(k + d) 

dU 2 (k ) 


= 2 p 2 0 (k)bg (k) + 2V (k) r 2 (k) > 0. 


(26) 


Noting the above equation, we know that the optimal control u (k) exists and can be determined by equating (25) to 
zero. Substituting (19) into (25) and let it equate zero we have 

P(k,q- l )A-\k,q^)B(k,q-'Mk)+[L(k,q-') 

+ P(k,q~ l )A~ 1 (k,q~ 1 )G(k,q~ 1 )]w(k) 

-Q(k,q-')z(k)+ r »° ( ' k)Vik) R(k, q~ l )u(k) = 0. 

Pno(k)b 0 (k ) 

Solving for u(k) we have 

[P{k,q-')A-\k,q-')B(k,q-')+ r "> {k) V {k > Rjk , q~' )}u{k) 

P-oWoW 

= Q(k,q~ l )z(k)-[L(k,q- 1 ) !JS) 

+ P(k, q~' )A _1 ( k , q~' )G{k, q-' )] w(k). 
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Left dividing (28) using P(k,q ')A '(k,q 1 ) we have the optimal control law as follows. 

[B(k,q ') + A{k,q ')P-'(k,q- ] ) R(k , q '^(k)) 

Pno (k)b 0 (k ) 

= A(k , q 1 )P~' (k, q' )Q(k,q 1 )z(k) (29) 

- [A(k, q 1 )P ik, q 1 )L(£, < 7 '') +G(£, 4 ~‘ )MD 

Where the noise w(k) can be estimated based on (1) using (21). Replacing w(kj using its estimate we have the GMVC 
controller (22). Comparing (1) with (21) we have 


C(k,q l )w(k + d) = 0 


(30) 


Where 


w(k + d) = w(k + d) -w(k + d) (31) 

Is the estimation error that will always decay exponentially to zero regardless initial conditions because of 
exponential stability of C~ l ik,q~ l ). From (22) and (31) we have 


[Bik, q~ l )+A(k, q-')p- l ik, q-') n0 )’\ R(k, q~ l )]uik )) 

P„ 0 ik)b 0 ik) 

= A(k,q- l )p-\k,q- 1 )Q{k,q- 1 )z(k ) 

-[A{k,q- l )p-\k,q- l )Uk,q~ l ) 

+Gik, q~ l )][w(k) - wik) + wik)\ 

= A(k,q- l )p-\k,q- l )Q(k,q- l )z(k ) 

+ [ Aik, q~ l )P~ 1 ik, q- 1 )Uk, q~ l )+ Gik, q- 1 )] [wik) -wik)]. 

Noting (1), (31) and (32) we have the following closed-loop equation for the LTV GMVC. 


Cik — d,q ’) 

0 

0 

wik) 

-Eik.q-') 

Tik.q-') 

0 

uik) 

0 

- Bik - d,q~')q~ d 

Aik - d,q~') 

y(k) 


0 

- Eik,q _1 ) 
Cik - d,q ~ l ) 


0 

A(k, q~')P~' ik,q~' )Q(k, q ~' ) 
0 


wik) 

z{k) 


(33) 


Where 

Tik,q-') = Bik,q “‘) 

+ A{k,q x )P l (k,q~ l ) 


ho-U r ll0 (k)V(k) 


Rik,q~ l ) 


P„ o(k)b 0 ik) 

Eik, q -') = Aik, q )P 1 (k, q^)Lik, q ~ l ) + Gik, q~' ). 


(34) 
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The closed-loop equation of the LTV GMVC shows that the closed-loop stability is determined by the following 
LTV matrix ARO. 


C{k—d,q~ l ) 0 0 

= - E(k,q "’) T(k,q ~ l ) 0 

0 - B(k — d ,q~ 1 )q~ d A(k-d,q- 1 ) 


(35) 


Noting that M(k,q~') is lower triangular and both A~\k,q~ l ) and C ' (k.q 1 ) are exponentially stable we know 
that the closed-loop system is exponentially stable if and only if T~ x (k,q~ x ) is exponentially stable. Equation (34) shows that 
exponential stability of T~ l (k,q~ l ) can be modified by choosing the weighting V(k) and the LTV filters P(k,q ' ) and 

R(.k,q~ l ). 


SIMULATION 


We consider the first order SISO LTV CARMA model 

y(k+2)+a(k)y(k+l )=u(k)+b(k)u(k-l )+w(k+2)+c(k)w(k+l ) (36) 

Where w(k) is an independent, stationary and Gaussian white noise with zero mean and unit variance. The plant 
parameters are 


f 0.4(1 -0.8c"*) 20/ <k< 2(i + 1) 

a(k) = < 

[—0.4(1 — 0.8e -t ) 20(i - 1) < ft < 20/ 
b(k) = 1 .3(1 + sin(0.4;zfc)) 


c(k) = < 


0 . 6 - 


- 0.6 


k + 1 
k 


k + 1 


10/ <k< 10(/ + 1) 
10(/-1) < k < 10/ 


(37) 


Where i=0 , 1, 2, We choose V(k)=l and 

P(k, q _1 ) = Q(k, q ~' ) = (1 + 0.2^ ‘)(1 + 0. l^ 1 ) -1 
R(k, cf') = (1 - 0.3<7 1 )(1 + 0.5ry~' ) _1 . 


Figure 1 shows the time-varying plant parameters. Figure 2 illustrates the plant output and the reference showing the 
LTV GMVC is able to stabilize the system and drive the output to follow the reference. Figure 3 shows the control variable 
u(k). 
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CONCLUSIONS 

An LTV GMVC has been developed for SISO LTV stochastic systems. It employs LTV filters in the form of moving 
average and autoregressive operators for the plant output and input providing a unified approach of LTV GMVCs for SISO 
CARMA models. It can be applied to a wide class of LTV systems. 
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